Restoration of broken symmetries in Self-Consistent RPA 
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It is shown that the Self-Consistent RPA (SCRPA) approach allows in a very natural way to 
restore symmetries, spontaneously broken on the mean field level. This is achieved via the intro- 
duction of a second Lagrange multiplier which constrains the variance of the symmetry operator to 
zero. This important feature of SCRPA, here pointed out for the first time, is illustrated employing 
a simplified model of the nuclear superfluidity. 



PACS numbers: 21.60.Jz 21.60.-n 21.60.Fw 71.10.-w 

Extensions of RPA theory, based on the Equation of Motion (EOM) approach have by now a quite long history. 
They, to a great deal, have been developed in nuclear physics. It started out with the work of Hara who included 
ground state correlations in the Fermion occupation numbers Q| . More systematic was the consequent work by Rowe 
and coworkers (see the review by D. J. Rowe Q). The same theory was developed using the Green's Function method 
by one of the present authors in |3j. Independently the method was proposed by R. Zimmermann and G. Ropke plus 
coworkers using a graphical construction Q . These authors called their method Cluster-Hartree-Fock (CHF) but it 
is equivalent to Self-Consistent RPA (SCRPA). The latter approach has recently been further developed by Dukelsky 
and Schuck in a series of papers jE] However, also other authors contributed actively to the subject 0. A number 
of remarkable results have been obtained with SCRPA in non trivial models where comparison with exact solutions 
was possible HJl(}. F° r instance for the exactly solvable many level pairing model of Richardson Jll[ SCRPA provides 
very accurate results for the ground state and the low lying part of the spectrum [0 . 

One important problem which remained little explored so far is the question how to deal with situations when 
there appears a continuously broken symmetry on the mean field level like nuclear deformation or superfluidity. We 
know that the standard RPA restores the broken symmetry in conjunction with a self-consistent mean-field single 
particle basis [Q. Complying with the Goldstone theorem is very important, since this assures simultaneously that 
conservation laws and sum rules are fulfilled. Unfortunately, it is a notoriously difficult problem that these properties 
are not easy to maintain when going beyond the standard RPA approach. SCRPA can be derived variationally and it 
couples in a non linear way the ordinary mean field to the quantal fluctuations J^|^] . If a symmetry is spontaneously 
broken it is customary in the Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB) approaches to fix at least the 
mean value of the symmetry operator with a Lagrange multiplier. There have been some attempts, with not very 
convincing outcome, to fix also its second moment (see for example the method proposed by Lipkin and Nogami in 
the case of pairing p3| ). HF and HFB approaches lead to one body mean field equations and it does not seem very 
natural to fix at the same time the mean value of a one body and a two body operator, that is e.g. the second moment 
of the symmetry operator. On the other hand in SCRPA which couples one body and two body mean field equations 
(see below) the requirement of simultaneously imposing the correct first and second moment is quite natural. This 
idea has to our knowledge not been explored so far and we present it here for the first time for a very simple model 
case of nuclear superfluidity. We will develope the general theory in a future work but in this short letter we think it 
appropriate to work out the theory in the very transparent case of the seniority model (l^j where we indeed will show 
that SCQRPA, as standard QRPA, completely restores the symmetry. However, contrary to standard QRPA, in the 
case of SCQRPA, because it is variational, this is equivalent to solving the model exactly. 
The Hamiltonian of the seniority model is given by : 



where G is the pairing strength, 2Q is the degeneracy of the level and the quasi-spin raising and lowering operators 
are given by : 



with m the quantum numbers of the magnetic substates. Since we will work with the quasiparticle RPA (QRPA), we 
will use instead of (pi) the following constrained Hamiltonian: 
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H' = H-mN- fi 2 N 2 (3) 

where N — X)m>o( a m a " l ~'~ a -ni a - m ) ^ s ^ e P ar ticle number operator, fix and /12 are two Lagrange multipliers allowing 
to fix the mean values N — (N) and iV 2 = (N 2 ) as we indicated above. We will transform H' to usual Bogoliubov 
quasiparticles Q : 

din = Ua^ + TO_ m , U 2 + V 2 = 1 (4) 

and define new quasi-spin operators as : 

P* = -J= 5>J„a! m , P=(P^ (5) 

^ m>0 

where the quasiparticle number operator is given by : 

= E + al m a_ m ). (6) 
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We will not repeat here how the Hamiltonian H' is expressed in the quasiparticle operators, since this is given in the 
text books [jl2). However, we briefly want to outline the main steps of the SCQRPA approach. The RPA excited 
states are, as usual, obtained as : 

\v) = C$\KPA) (7) 

where \RPA) is the correlated RPA ground state defined via the vacuum condition : 

Q\RPA) = 0, (8) 

and the QRPA excitation operator is given by ||] : 

Q^=xP^-yP, x 2 -y 2 = l. (9) 

It is directly checked that relation ([}]) can be inverted, that is : 

P f = .xQt + yQ (10) 

and therefore with condition (||) a certain number of expectation values are readily expressed by the RPA amplitudes 
x and y : 

(PPt) = x 2 {1 _ <M), 
(Ptpt) = (PP) = xy(l - M) (11) 

where (...) = {RPA\...\RPA)/{RPA]RPA). With (g) and (|) it is standard procedure using e.g. the EOM approach 
to derive the RPA secular equation |2] : 



or explicitly : 
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It is important to note that eqs ( |l2f]lq ) can also be derived from a variational principle using the following functional 
flit , which is a special combination of the energy weighted sum rule JlJ] : 

[ ,yJ EJHQ t |o)| 2 -E,'I^IWI 2 [ ] 

with i?^ and \v) eigenvalues and eigenstates of H' . 

Since (^) can be considered as a Bogoliubov transformation of Fermion pair operators we can interpret the SCQRPA 
equations as the corresponding mean field equations. In other words the SCQRPA equations are the mean field 
equations for a gas of quantal pair fluctuations. The standard mean field equations of the one body type are, as 
usual, derived in minimizing the constrained ground state energy (RPA\H'\RPA) = E' with respect to the mean 
field amplitudes u, v : 

- B ± = {[H',P\) = 0. (17) 

This last equation (which can be verified in evaluating l.h.s and r.h.s separately) is very interesting, since it is equivalent 
to ([H',Q]) = i-§t(Q(t)) = which should be fulfilled at equilibrium and consequently constitutes a further relation 
of the EOM approach, first put forward in H. With the normalization of the amplitudes u, v (see (^)) one can readily 
rewrite ( |l7| ) in the more common form : 
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where, 

h = -{G + A^i 2 Sl)v 2 - fix, (19a) 

f ft \ _ <£i> l 

A = (GSl + 4 M2 )m, - (G - 4(i 2 )uvl 2(xy + y 2 ) + X ql ^ L (19b) 



= \/h 2 + A 2 . (19c) 



We see that once we replace in (19) \RPA) by the uncorrelated HFB vacuum with 

V = (n q ) = (n 2 q ) = and n 2 = 



we obtain from ( |18| ) the standard HFB equations. However, in keeping with the correlated vacuum the Fermionic 
mean field equations ([l^) get coupled to the mean field equations (|l^) for the quantal pair fluctuations which is a 
very gratifying aspect of SCRPA theory. 



In addition to the two mean field equations (13) and ( |18[ ) we have two further equations which, in principle, allow us 
to find the Lagrange multipliers fix, \i 2 (see, however, below) : 

N = (N) = (u 2 ~ v 2 )(h q ) + 2Slv 2 , (20a) 

N 2 = (N 2 ) = (u 2 ~v 2 ) 2 {h 2 q )+8Slu 2 v 2 (l- { -^)(xy + y 2 )+4S}v 2 (u 2 + Slv 2 ) 

+ Av 2 (SI(u 2 -v 2 )-u 2 ){n q ) " (20b) 



We again see that eqs (|20|) reduce to the standard expressions, once, as in the H FB a pproximation, we pose y = 
(h q ) = (ft 2 ) = 0. In the particular case of the seniority model the number equation (20a) in the HFB approximation 
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determines the amplitudes u,v and then no freedom is left to impose AN = y (N 2 ) — (N) = 0. However, in the 
present more general approach there is more freedom and we will be able to satisfy the relation AiV = 0. Before 
getting to this point, we first have to show how the so far open quantities (n q ) and (h 2 ) can be expressed via our four 
unknowns v, y and u 2 (the two other amplitudes u and x being determined via the corresponding normalization 
conditions in (||) and The determination of those two quantities is one of the difficulties in the SC(Q)RPA 

approach [^|^]. However, this problem has found an elegant solution in the early work of ]l5[ |. There is no place 
here to go into the details of the derivation. Let us simply state that applied to the present model one is able to 
expand n q = ^ i=0 Ci(P^) l {P) 1 and also n 2 = J2i=obi(P^Y(PY where the coefficients Ci,bi have definite expressions. 
Explicitly this yields [|l6| : 

- n p,f + mhr> pi ' 2p2 + n» ( n - h - ^ * p * + - <21a) 

^-^m-^^w^rr^ 



We see that between (21a) and (21b) exists, at any step of the expansion, the following relation : 



-(pip + P P i) + (^-—f = -(- + !). (22) 



It should be noted that eq (22) is just the Casimir relation for the quasi-spins of the seniority model. Using (10), (|8|) 
we obtain from ((2^), ( p2|) t wo eqs for (n q ) and (n 2 ) which can be solved in terms of v and y amplitudes and therefore 
the two mean field eqs (|13|) and (113) together with (|20aD a nd (2 0b) represent four coupled nonlinear equations for our 



unknowns v, y, \i\ and /i 2 . Usually the number equations (20a) and (20b) are to be used for the determination of the 



chemical potential \i\ and the second Lagrange multiplier /i 2 and the mean field equations for the amplitudes u, v and 
x,y. In the present case it is, however, more convenient to invert the role of mean field and number equations, since 
eqs ( pp[ ) do not depend on the Lagrange multipliers and therefore readily allow to determine v 2 and y 2 as a function 
of the particle number N. Inversely the two mean field eqs (|l3|), ( |Tg| ) are linear in /ii and /i 2 and for instance it is 
seen that ( |l3| ) directly yields : 

»2 = ~ (23) 

independent of the particle number N. Considering, the well known exact expression for the ground state energy of 
the model @: 

E = -^(n + l)N+^N 2 (24) 

we see from ^ 2 = f |p- t hat (|23|) gives the exact value for the second Lagrange multiplier ^ 2 ■ For the chemical potential 
/ii we obtain from ( jl8| ) 

r 1 if, \ — ~\ n 

Hi = (G- 4 M2 ) (n l)v 2 + (1 - 2v 2 )(xy + y 2 + - S q, _ (fl * ) - -Q - 2/i 2 . (25) 

With relation (^3|) this gives pb\ — —-^(O + l) which again is the exact value. Furthermore, with (|2^) we have from 
(11), (0) that A = B = and therefore the RPA eigenvalue E = 0. This means that, as in standard QRPA, SCQRPA 
yields a Goldstone mode at zero energy which in nuclear physics usually is called the spurious mode jl2| . This feature 
is very rewarding, since it signifies that the particle number symmetry is exactly restored. 

It is well known that restoration of good particle number implies in this very simple model case that the model is 
solved exactly J^]. We have already seen that one obtains the exact values for /ii and /i 2 . We now will show that one 
also obtains the exact value for the ground state energy (and therefore for the whole band of ground state energies) . 
This goes as follows. For the expectation value of H of eq (|l|) in the RPA ground state, using (||), ( |Io[ ) and the 
quasiparticle representation for H, we can write 



E = (H[ 

= -|(n + !){(! - 2v 2 )(h q ) + 2Uv 2 } + 2v 2 ) 2 (h 2 q ) 



4 



+ [Av 2 (n(l - 2v 2 ) -l + v 2 )- 8(1 - v 2 )v 2 (xy + y 2 )](h q ) 

+ 8fi(l - v 2 )v 2 (xy + y 2 ) + 4f7w 2 (l - v 2 + ilv 2 )]. (26) 



In this expression we have used the relation 4j/ 2 (f7 — (rig)) = 2(f2+ !)(«■<?) — (n 2 ) which follows from (123). Using the 



expression for N and N of (20a) and (20b) once more, we see that the exact expression 124) is recovered. 



Before concluding several remarks are in order. One may have noticed that the model is solved exactly independent 
on how far we push the expansion (^T|). This is, of course, particular to the simplicity of the present model. As long 
as the theory is sufficiently general to fulfill the two number equations ( p0a| ) and (20b), the model will always be 
solved exactly in the present approach. One also notes that we have developed our theory without ever using the RPA 
ground state wave-function explicitly. Equation (^) was used only implicitly. In the present case the vacuum condition 
can be solved and an explicit expression for \RPA), quite analogous to the ones found in other simplified models [^), 
can be given. We will not present this here, since it turned out in more realistic cases that the explicit construction 
of the RPA ground state is extremely difficult, if not impossible. We would like to point out, however, that also in 
HFB theory, in principle, the explicit knowledge of the wave function is never necessary and for all relevant questions 
it is sufficient to deal with the vacuum relation a\HFB) = 0. 

In conclusion we showed in this work how in the case of a symmetry broken situation the introduction of a second 
constraint which yields the correct value (zero) for the variance of the symmetry operator allows to restore exactly 
the symmetry, and thus the conservation laws, in the SCRPA approach. This was demonstrated using a simplified 
pairing model. It turns out that this solves the model exactly. Far from being over simplified for the present purpose 
we think, on the contrary, that the model perfectly illustrates the flexibility of the SCRPA to accommodate the ad- 
ditional constraint of putting the variance of the symmetry operator to zero. Indeed in this model, precisely because 
of its simplicity, there is within HFB no freedom for a constraint other than the single one for the mean value of the 
particle number. We think that the possibility of fixing simultaneously mean value and the variance of the symmetry 
operator is a very important extension and ingredient of SCRPA. This has been worked out for the first time in this 
present note. It implies that conservation laws and usual sum rules remain fulfilled in the same way as in the standard 
RPA. This is a very remarkable feature of SCRPA, since the conservation of symmetries is a well known problem 
when attempting to include correlations beyond standard RPA. The SCRPA approach, contrary to the standard RPA 
which is perturbative, can be derived variationally and is thus non perturbative. We intend to apply the present 
theory to more complicated situations and to give a more elaborate account of it in a future publication. 
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